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ABSTRACT 



We present an analysis of the broadband UV and optical properties of 2; < 3.4 
quasars matched in the Galaxy Evolution Explorer (GALEX) General Data Re- 
lease 1 (GRl) and the Sloan Digital Sky Survey Data Release 3 (SDSS DR3). Of 
the 6371 SDSS DR3 quasars covered by 204 GALEX GRl tiles, 5380 (84%) have 
near-UV detections, while 3034 (48%) have both near-UV and far-UV detections 
using a matching radius of 7". Most of the DR3 sample quasars are detected 
in the near-UV until z ~ 1.7, with the near-UV detection fraction dropping 
to ~ 50% hj z ~ 2. Statistical tests performed on the distributions of non- 
detections indicate that the optically selected quasars missed in the UV tend 
to be optically faint or at high redshift. The GALEX positions are shown to 
be consistent with the SDSS astrometry to within an rms scatter of 0.6-0.7" in 
each coordinate, and the empirically determined photometric errors from multi- 
epoch GALEX observations significantly exceed the Poissonian errors quoted in 
the GRl object catalogs. The UV-detected quasars are well separated from stars 
in UV/optical color-color space, with the UV/optical relative colors suggesting a 
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marginally detected population of reddened objects due to absorption along the 
line of sight or due to dust associated with the quasar. The resulting spectral 
energy distributions (SEDs) cover ~ 350 — 9000 A (rest frame), where the overall 
median SED peaks near the Ly a emission line, as found in other UV quasar 
studies. The large sample size allows us to construct median SEDs in small bins 
of redshift and luminosity, and we find the median SED becomes harder (bluer) 
at UV wavelengths for quasars with lower continuum luminosity. The detected 
UV/optical flux as a function of redshift is qualitatively consistent with attenu- 
ation by intervening Lyman-absorbing clouds. A table containing the 6371 DR3 
quasars covered by GALEX GRl and the results of the matching is available in 
the electronic edition of the journal, and a table containing the overall median 
SED is also available. 

Subject headings: catalogs — surveys — galaxies:active — quasars:general — 
ultraviolet igalaxies 



1. INTRODUCTION 



One of the primary goals of large-scale UV and optical sky surveys is to create large 
sample sizes, crucial for a statistically significant measurement and interpretation of the 
ubiquitous properties displayed by active galactic nuclei (AGN), such as broad emission 
lines and power-law continua. In quasars, a large fraction of the bolometric luminosity is 
emitted in a strong, broad feature that begins to dominate the spectral ener gy distribution 



(SED) at the bluest optical wavelengths (the 'big blue bump'; hereafter BBB; ISanders et al. 



19891 ) and appears to extend shortward to the current limits of UV satellite data (~ 100 A). 
Previous work that interprets the rise in the continuum as the thermal signat ure of a large 



and hot accretion disk surrounding a central supermassive black hole (e.g., IShieldd Il978 



Malkan fc Sargentlll982l ) has argued that the BB B should peak in the extreme ultraviolet 



(EUy), shortward of th e Lyman limit at 912 A ( Mathews fc Ferlandl 119871 : iHubeny et al. 



2000l : iFrank et al.ll2002l ). Superimposed on the BBB is yet another feature (the 'small blue 
bump'; hereafter SBE i) attributed to Balmer contiri uum and complex Fe II emission from 
the broad line region (jGrandilll98ll : IWills et al.lll985l ). The BBB and SBB are characteristic 
features of quasar emission; therefore UV and optical observations are a critical probe of the 
physics of the inner regions of quasars. Deep UV observations have the added capability 
of sampling the ionizing flux, which gives rise to the powerful emission lines that are a 
trademark of quasar spectra. 
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Zheng et al.l (119971 ) have used Hubble Space Telescope (HST) observations to construct 
composite rest-frame UV spectra from 101 AGN with z > 0.33 that display a break in the 
UV/optical power-law continuum near ~ 1000 A (specifically, the EUV continuum blueward 
of the break is soft er than the continuum redward of the break). This result was confirmed by 
Telfer et al.l (120021 ) for a larger sample of 184 AGN with redshifts in the range 0.33 < z < 3.6, 
who also demonstrated no evolution of the EUV spectral index a {fy oc z/") w ith redshift 



and o nly a moderate trend with luminosity for the radio-loud objects. However, IScott et al. 



(120041 ) concluded that the UV continuum from their composite of 85 AGN spectra with 
z < 0.67 taken from Far Ultraviolet Spectroscopic Explorer {FUSE) archival data, can be 
characterized by a single power law. Given that the FUSE sample spectra have enhanced 
emission lines relative to the HST sample and have a median luminosity that is nearly an 
order of magnitude lower, this suggests that the strength of any continuum break may be 
luminosity dependent. The results obtained from the redshift and luminosity distributions 
of the FUSE and HST AGN samples also suggest possible evolution of the UV SED with 
redshift. 

High S/N rest-frame UV/optical composites have also been constructed from large 
samples in ground based optical surveys, such as the com posite of 688 optic ally-selected 
0.2 < z < 3.0 quasars from the Large Bright Quasar Survey (IFrancis et al.lll99ll ). 657 radio- 
selected 0.23 < z < 3.40 q uasars (^with a larger fr action of low-redshift objects) from the 
First Bright Quasar Survey ( IBrotherton et al.ll200ll ). and for 2200 optica lly-selected quasars 



spann ing 0.04 < z < 4.80 from the Sloan Digital Sky Survey (SDSS; iVanden Berk et al. 
200ll ). These studies have provided great insight into the detailed properties of quasars 



longward of Ly a, but optical spectra alone can only probe the regions shortward of Ly a 
A1216 at higher redshifts {z > 2), where the density of Lyman-absorbing clouds is high and 
significantly limits both the number of detections and the S/N. 

We investigate the results reported for FUSE and HST UV composite SEDs by ex- 
panding the sample size by more than an order of magnitude — covering the full range of 
luminosity and redshift spanned by both studies. We also perf orm an empirical analysis of 
the Galaxy Evolution Explorer ( G ALEX ■ AMoirissey et al.ll2005l ) astrometry and photometry 
that is crucial for any assessment of the reliability of GALEX positions and interpretations 
of the resulting UV/optical colors and SEDs. In §2], we describe the data sets used to find 
UV/optical candidates taken fror n the GALEX G eneral Data Release and the 



Sloan Digital Sky Survey (SDSS; lYork et al.l l2000l ) Data Release 3 (DR3; lAbazajian et al. 



20051 ). as well as the description of the construction of a clean and statistically significant 
photometric sample of quasars. An analysis of the UV/optical colors of the photometric 
sample of GALEX /SDSS matches, as well as the resulting SEDs, are given in ^ and ^ 
summarizes our results obtained from the clean sample and the interpretation of the mean 
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quasar SED luminosity and redshift dependence. Throughout this paper we assume a A- 
dominated flat cosmology with Q\ = 0.7, 0^ = - 3, an d Hq = 70 km s"^ Mpc~^, consistent 
with the latest WMAP results of ISpergel et al.l (120061 ) . In our study, none of the objects 
in our UV/optical sample have been corrected for photometric contributions from the host 
galaxy. Unless otherwi se indicated, we corr ect all photometry for Galactic extinction using 
the extinction maps of ISchlegel et al.l (Il998l ). 



2. THE DATA SET 

The inclusion of UV observations allows for (1) diagnostics that can probe the evolution 
of the UV/optical SED with redshift, and (2) observations of the rest-frame UV emissio n from 



lower redshift quasars that are not heavily affected by the Ly a forest of absorption (ILynds 



19711 ). which one cannot accomplish with optical surveys alone. The accu rate photometry 



0.02 mag) and size (46,420 quasars) of the SDSS DR3 Quasar Catalog (ISchneider et al. 



20051 ) provides a unique context for the interpretation of UV observations. 



2.1. GALEX 



The GALEX General Data Release 1 fGRl: iMorrissev et al.ll2005l : Martin et all 120051 1 
includes observations covering nearly 8000 deg^ of sky in four imaging surveys and two grism 
surveys as of December 2004. GALEX is a NASA Explorer Mission that is performing an 
All-Sky Imaging Survey (AIS) to itlab ~ 20.5, a Medium Imaging Survey (MIS) covering 
1000 deg^ to niAB ~ 23, a Deep Imaging Survey (DIS) covering 100 deg^ to niAB ~ 25, and 
a Nearby Galaxy Survey (NGS) with a surface brightness limit of niAB ~ 27.5 arcsec"^. The 
observations are obtained in two UV bands — the far-ultraviolet (the FUV band) covering 
~ 1350 — 1750 A and the near-ultraviolet (the NUV band) covering ~ 1750 - 2750 A. The 
AB-magnitudes corresponding to the FUV and NUV bands will hereafter be denoted by / 
and n, respectively. The 1?25 field-of-view was designed to attain 4"-6" resolution; currently 
the pr ocessing pipeline is achieving a FWHM of 6'.'5 in the FUV band and 7'.'2 in the NUV 
band (IMartin et al.ll2005l ). The large sample of objects provided by GALEX, along with 
deep surveys in the optical and the near-IR, provides multiwavelength coverage for a large 
number of galaxies, constraining star-formation models and their evolution with cosmic time. 
Lyman absorption causes the flux of many optically identified sources to quickly diminish 
with redshift at shorter wavelengths; however, the extended coverage beyond the rest-frame 



^GRl documentation can be accessed at |http: / /galex. stsci.edu/ GRl / 
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Lyman limit at 912 A for higher redshift objects has allowed for a number of UV counterparts 
to optical SDSS sources to be identified. For many objects, GALEX also allows a glimpse of 
the wavelength regime shortward of Ly a, and has permitted observations as blue as ~ 430 
A in the rest frame {z < 2), beyond which the Lyman absorption becomes very strong. 



2.2. SDSS 



The SDSS will map ~ 10^ deg^ of the sky (mainly concentrated at high Galactic lati- 
tudes) and obta in spectra for ~ 10^ galaxies and ~ 10^ qu asars in five opti cal photometric 



bands (n.q.r.z.z.' iFukugita et al.lll996l ) to a depth of r ~ 22 (lYork et al.ll2000[) — with photo- 



metric errors of less than 0.02 magnitudes Jlvezic et al. 2004 ). and an astrometric accuracy 
better than O^^l for r < 20 5 JPier et all boosh. The SDSS obtains im aging with a CCD 
mosaic camera (iGunn et al.lll998l ) coupled to a dedic ated 2.5m telescope (|Gunn et al.ll2006r) 
at the Apache Point Observatory. The SDSS DR3 ( Abazajian et al.ll2005l ) contains more 
than 10^ cataloged objects covering over 5000 deg^ of sky, and the Q uasar Catalog include s 
observations of 46,420 type 1 quasars with > —22 and z < 5.41 (ISchneider et al.ll2005l ). 
In common with the GALEX phot ometric system, the u, q, r, i, and z magnitudes are 



normalized to lie on the AB-system (jOke &: Gunnlll983l : ISmith et al.ll2002l ). 



The SDSS quasar survey is largely motivated by two science goals: (1) to address 
the evolution of the quasar luminosity function with time and (2) to exami ne the redshift 



depen dence of quasar clustering. The SDSS quasar target selection algorithm (IRichards et al. 



2OO2I ) uses point-spread function magnitudes corrected for Galactic ex tinction; it iriyolves 
targeting those objects that are unresolved in the FIRST radio survey (IBecker et al.lll995l ) 
and occupy distinct regions of mult i- dimensional color space outside of the stellar locus and 
the regions occupied by nonactive galaxies — optimized for the low- and high-redshift objects 
and to detect quasars for < z < Q. For quasars with i < 19.1, ~ 95% are found by the 
SDSS color-selection technique (IVanden Berk et al.ll2005l ). The DR3 Quasar Catalog (46,420 
objects) consists of images taken between 1998 September 19 and 2003 May 1, as well as 
spectra obtained between 2000 March 5 and 2003 July 6. The bright limit of the target 
selection algorithm is i = 15.0. 



2.3. Cross-Matching and Astrometry 



Previous work finding matches between GALEX and the SDSS has been implemented 
using GALEX Early Release Observations (ERO) for fields not restricted to quasars (e.g.. 
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Agiieros et al.ll2005l ). for internal data rele ases fe.g.. iBianchi et al.ll2005l : ISeibert et al.ll2005l ) 



and for small regions covered by Spitzer (IRichards et al.ll2006l ). We searched for UV coun- 



terparts to the 6371 SDSS DR3 quasars covered by 204 GALEX tiles in the GALEX GRl 
merged object catalogs (hereafter referred to as the 'quasar sample') — using an initial search 
radius 9 = 7" (roughly 1" larger than the FWHM of the angular resolution limit in the NUV 



band; for a more complete description see iMorrissey et al.ll2005l ). For a mean SDSS quasar 



separation of ~ 20' on the sky, the astrometry of GALEX (see below) and the SDSS is more 
than sufficient such that the typical astrometric errors for a given source location will be sig- 
nificantly smaller than the average quasar separation — this simplifies the matching routine 
significantly, since the probability of more than one quasar lying within our search radius is 
extremely small. In the case of more than one GALEX UV source coincident with an SDSS 
DR3 quasar to within 7", the closest one was selected. When overlapping GALEX GRl 
tiles each covered the same SDSS quasar, we have selected the observation from the tile with 
the highest effective exposure time (tefr) in the NUV band (and, in general, correspondingly 
smaller photometric uncertainties for the same objects due to higher S/N). 

After matching the DR3 quasar sample with the GRl catalogs, 5380 quasars (84.4%) 
are detected in the NUV, while 3034 (47.6%) are detected in both the NUV and the FUV. 
We also note that 38 quasars (0.6%) were detected only in the FUV (see discussion below). 
There were 953 (14.9%) quasars with no GRl counterparts; a statistical analysis of the DR3 
quasars with and without a UV detection is given in §3.11 The mean separation between a 
DR3 quasar and the nearest GALEX UV source was ~ 1", with ^ 25% of the quasars having 
more than one UV source within 7". If a specific tile has been observed multiple times, then 
we use the coadded result; however, the GALEX source extractor can and does identify UV 
sources unique to each epoch, so it is possible that an appreciable fraction of the quasars 
with multiple UV sources within 7" are simply due to this phenomenon. For the quasars 
with multiple associated UV sources, the mean nearest separation is 0'.'7, while the mean 
separation of the next nearest UV source is 1'.'4. 

Our matches with only an FUV detection were rare — with 18 sources appearing on 
MIS tiles, 19 on AIS tiles, and a single source on a NGS tile. Diffuse reflections from bright 
stars located immediately outside the GALEX detector fleld of view are not flagged as image 
artifacts in the merged object catalogs, and can artiflcially increase the number of false point 
source extractions for a given fleld. In addition, the source extractor of the GALEX image 
processing pipeline will often fragment resolved objects into individual UV sources, thereby 
contaminating search results with multiple matches to the same optical For 
the remainder of our analysis, we exclude the FUV-only matches to the quasar sample, 
and given the small fraction of the total sample with FUV only detections, we expect the 
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remainder of our matches to be negligibly contaminated by false detections. 

Figure [T] illustrates the coordinate offsets between SDSS and GALEX positions for all of 
the matches in both UV bands. These results show that ~ 95% of our matches lie within a 
matching radius 6 = 1.7", while ^ 99% are included if = 2.6". The median offsets in right 
ascension and declination are -0"03 and -(y.'12 respectively. Monte Carlo simulations were 
used to generate predicted distributions of coordinate offsets between GALEX and SDSS po- 
sitions, using a set 10^ objects and assuming an rms ctsdss = O'-'l ioi each SDSS coordinate. 
A comparison of the simulation results with our observed distribution of positional offsets 
shows that our statistics are consistent with acALEX = 0.6-0.7" for each GALEX coordinate 
(the O'.'l rms for the majority of SDSS objects contributes negligibly when added in quadra- 
ture). Our sample is not significantly contaminated by false matches (there is no large excess 
of objects at larger separations); although, as one might expect, there are noticably more 
high-cr matches than predicted by our simulations, primarily because our Monte Carlo runs 
do not take into account possible systematic errors affecting the GALEX astrometry (e.g., 
the rate of false matches introduced by multiple source extractions of resolved sources or im- 
age artifacts). Our assessments of the G^lLi^X astrometry by making use of SDSS positions 
are fully consistent with the astrometry derived from ground-based calibrations conducted 
by the GALEX collaboratioi]f|. 

For a completely random distribution of GALEX UV sources in the deepest (DIS) 
fields, the probability of any (star or galaxy) UV source coincident with an SDSS quasar 
for e = 7" is ^ 10%, while it is ^ 1% for 6 = 2!!Q. Given that the majority (> 99%) of 
the SDSS quasars are covered by shallower tiles, we use 1% as a liberal upper bound to 
our contamination rate by non-quasar counterpart UV detections (either by stars or image 
artifacts). Accordingly, for the remainder of our analysis and our discussion we will refer to 
the 'detections' as those matches with GALEX /SDSS separations less than 2'.'6 (5165 NUV 
detections and 3004 NUV+FUV detections). 



2.4. GALEX Photometry 

In contrast to observations of relatively bright objects in the local Universe (e.g., nearby 
resolved galaxies, stars, etc.) where the photometric errors are usually small, photometric 
uncertainties of relatively faint objects such as quasars in flux limited surveys are no longer 



^For additional GRl documentation 
[http://wwwTgalex.caltech.edu/DATA/[ 



and 



description of data caveats 



see 



The GvlLiJX collaboration website is littp://www. galex.caltech.edu/ 
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negligible and can largely determine the confidence of the interpretation of the overall photo- 
metric properties of a sample. UV observations are more sensitive to line of sight extinction 
than optical observations; thus quasars (being fainter and at higher redshifts) are not as read- 
ily detected in the UV, despite emitting a significant fraction of their bolometric luminosity 
in the UV region. For the ones that are found in the UV, the photometric uncertainties 
can become comparable to the detected flux itself, and one must therefore know the distri- 
bution of errors well before forming conclusions (e.g., the resulting UV colors and SEDs of 
optically-selected quasars) in the context of the rest-frame emission. 

The GRl UV source extractor of the GALEX image processing pipeline quotes photo- 
metric uncertainties for each object by assuming the observations are Poisson noise limited 
(based on the number of detected photons). Therefore, the magnitude errors given in the 
GRl object tables do not take into account additional sources of noise — including unknown 
variances of the detector background level and flat-field maps, or any other systematic errors 
present in the data. If one assumes the bulk of the stars observed in a GALEX field are 
non-variable, then fields with multi-epoch observations are a useful tool for analyzing the 
repeatability of the GALEX UV photometry for the same objects, and one can use large 
numbers of objects to empirically estimate the true photometric uncertainties. In §31 we 
incorporate the results of using this method (described below) into our analysis of the UV 
colors of the quasars in our sample. 

Using multi-epoch GALEX tiles from the Spitzer First Look Survey area with similar 
exposure times taken from the Deep Imaging Survey, we matched 7090 NUV detections using 
a search radius of 7" (~ FWHM in the NUV band) and restricted the matching out to 95% 
of the tile radius, thus minimizing the contamination by image artifacts and multiple source 
extractions of diffuse out-of-field bright stars that can dominate the sources near the edge of 
the fields. After removing matches with magnitude differences > 3cr from the median of each 
magnitude bin, 6965 (98.2%) objects remain. Figure [2] shows the NUV magnitude differences 
(ANUV) as a function of NUV magnitude, with the red curves showing the errors given in 
the GRl tables after adding the errors in quadrature for each epoch (and reflected about 
ANUV = 0). The green and blue curves are the rms and 68.3% confidence limits respectively, 
which were calculated in NUV magnitude bins. There is good agreement between the rms 
and confidence limits, indicating that the errors in the NUV photometry are approximately 
Gaussian distributed at the locations of each magnitude bin. A similar analysis was also 
performed for repeated FUV observations, as well as for both the NUV and FUV bands in 
the All-Sky Imaging Survey (e.g.. Figure [3]). 

Figure H] shows the least-squares fits to the binned NUV and FUV magnitude rms as 
a function of magnitude. Open and closed symbols are calculated from the DIS and AIS 
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results, respectively, as discussed above. Due to the mean exposure time for an AIS tile 
being roughly an order of magnitude lower than that for a DIS tile, the AIS empirical error 
best-fit curves lie above those for the DIS due to lower S/N. The qualitatively larger error 
and fewer bins in the FUV results from the lower density of FUV sources on the sky above 
the threshold flux and the lower sensitivity of the FUV band relative to the NUV band. 
Figure [5] shows the distribution of effective exposure times for all of the quasar matches 
in both UV bands. Given the bimodal distribution, we divide the sample at an effective 
NUV exposure time of log (tofr) =2.7 — for our empirical estimates of the true photometric 
uncertainty, we employ the AIS curves for detections with log (tofr) < 2.6 and the DIS curves 
for log (teff) > 2.6 (see Figure H]). If we define a and b to be the slope and y-intercept of the 
overplotted curves respectively, then the empirical photometric error (Am) is given by 

log Am = am + b, (1) 

where m is the quasar NUV/FUV magnitude. Table [1] gives our best-fit values for a and 
b, which we have used to determine the empirical photometric errors given in our table of 
matching data (§2.5). 

From the number of repeated NUV detections across multiple epochs for a DIS field, we 
estimate an ~ 8% probability that a quasar from the SDSS DR3 will have have a false match 
in the NUV band for a search radius of 7", while this drops to < 1% if the search radius is 
limited to 2.6" — this justifies the selection process discussed previously in §2.3[ In addition, 
we extend this contamination rate to the entire quasar sample, given that these calculations 
are based on the deepest GALEX fields with the highest density of UV sources. The results 
of the empirical analysis of the GALEX NUV and FUV photometry clearly show the true 
photometric error rises above the Poissonian errors quoted in the GRl tables (e.g.. Figure 
[2]), where for the NUV band quasar detections (median n ~ 19.8), the photometric errors 
(assuming the stars to be non-variable) are already in excess of ~ 5% — further justifying 
this analysis as necessary for constraining the behavior of the quasar UV colors with redshift, 
the relative colors (gS3]), and the UV/optical SEDs (gS3D- 



2.5. Data Table 

Our database (see Table [2]) of GALEX GRl matches to SDSS DR3 quasars contains all 
of the parameters listed in Table [3] (and described below) for each of the 6371 SDSS DR3 
quasars covered by GRl tiles. Each entry consists of selected parameters from both the 
GRl catalogs and the SDSS DR3 quasar catalog, as well as our measured quantities. Full 
versions of the GALEX GRl tables are publicly available at http://galex.stsci.edu/GRl/ 



- 10 - 



and the DR3 quasar catalog is available via the SDSS public website at http: / /www. sdss.org. 
Of the 6371 quasars in the table, 5380 (84.4%) are NUV detections and 3034 (47.6%) are 
NUV+FUV detections, found by matching coordinates and before any quality cuts are made 
(see §2.3p . The table also includes 38 (0.6%) FUV only detections — these objects are likely 
contamination by false matches (a mean separation of ~ 2.2") and we therefore do not refer 
to them as 'detections' for the purposes of our analysis. 

For the DR3 sample quasars without an NUV and/or FUV source within 7", the entries 
in columns 33 - 34 are the NUV teff and/or FUV t^s taken from the tile that covered the 
quasar — the GRl tables do not give upper limits in magnitudes. In principle, the relative tes 
would be the only parameter necessary to calculate upper limits for a given tile, but the non- 
uniform response of the GALEX detectors as well as changing UV background levels would 
necessitate calculating upper limits on a tile-by-tile basis. Future GALEX data releases may 
enable a more straightforward calculation of upper limits for each optically-selected quasar 
without a UV detection. 

Description of the table parameters: 

1. The SDSS DR3 object designation from the quasar catalog in the format SDSSJhh- 
mmss.ss-|-ddmmss.s, where 'SDSSJ' is omitted. 

2. The GALEX GRl object designation from the merged object catalog. The same ob- 
ject identified on multiple GALEX tiles can have more than one object id. Thus, one should 
always use object coordinates and not object designations to identify individual objects as 
such. DR3 quasars without a UV detection (see §2.3p will have a '-999' entry. 

3-4. The object coordinates (right ascension and declination; J2000.0) from the SDSS 
DR3 quasar catalog, given in decimal degrees. For a significant majority of the objects, the 
SDSS astromet ric accuracy of e ach coordinate is O'.'l, with the largest errors expected not to 



exceed 0'.'2 fsee lPier et al.ll2003[ ) 



5-6. The object coordinates (right ascension and declination; J2000.0) from the GALEX 
GRl merged object catalogs, given in decimal degrees. Our empirical analysis of §53] demon- 
strates the derived rms of ~ 0'.'6 - 0'.'7 for each GALEX coordinate. DR3 quasars without 
any UV detection will have a '-999' entry for both coordinates. 

7. The coordinate offsets between the GALEX and SDSS positions {GALEX - SDSS), 
given in arcseconds. Objects with separations > 7" or non-detections have a '-999' entry 

8. The quasar redshifts from the DR3. For any given object, the estim ated redshift 



uncertainty is < 0.01, except for Broad Absorption Line (BAL) quasars (e.g.. iTrump et al. 



2006|). 
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9-18. The DR3 PSF magnitudes and errors for all of the quasars (no Galactic extinction 
corrections applied), given in AB-magnitudes from the five SDSS photometric {u, g, r, i, 
and z) bands. The magr iitudes are those deriv ed from BEST photometry (or TARGET if 



BEST is unavailable; see ISchneider et al.l (120051 ) and references therein) 



19-24. The GALEX PSF magnitudes and errors (both as quoted in the GRl object 
catalogs and from our empirical analysis of the photometry), in units of AB-magnitude from 
both UV photometric bands. Th e effective wavelengths of the FUV and NUV channels are 



1528 A and 2271 A, respectively (IMorrissev et all 120051 1. 



25. The Galactic extinction in the SDSS-'u band {Au) derived from lSchlegel et al.l (119981 ). 
The extinction in the remaining SDSS bands can be expressed as Ag, Ar, Ai, and A^ with 
values corresponding to 0.736, 0.534, 0.405, and 0.287 times Au- For the GA LEX f and n 



bands, the corrections Af and An are 1.584 and 1.684 times A^, respectively (IWyder et al. 



20051 ) . All correction factors assume an Ry = 3.1 absorbing medium. 



26. Log Nh represents the logarithm of the neutral Hydrogen column d ensity along our 



Galac tic line of sight to each quasar, interpolated via 21 cm data (also see ISchneider et al. 



(120051 ) and references therein). 



27. The absolute z-band magnitude — calculated after correcting for Galactic extinction 
and assuming a power-law continuum index of a = —0.5 {f^, oc z/"), and a A-dominated flat 
cosmology with = 0.7, flm = 0.3, and Hq = 70 km s~^ Mpc^^ (same for this paper). 

28. The logarithmic monochromatic quasar luminosity (z/L,^) computed at 2200 A (rest- 
frame). A '-999' entry is given if there is no wavelength coverage. 

29. The same as column 28, but for 5100 A — only computed for the quasars where 
5100 A has not been redshifted out of the SDSS-^; band. A '-999' entry is given if there is 
no wavelength coverage. 

30-31. The A{g — i) and A(/ — n) relative colors ( §3.3p . defined as the g — i and f — n 
colors minus the corresponding median colors at the redshift of the quasar, given in units of 
AB-magnitude . 

32. Object distance from GALEX field center in decimal degrees. Note that for the 
case of objects with multi-epoch observations across multiple GALEX fields, the entry from 
the tile with the highest effective NUV exposure time ('NUV_WEIGHT' in the GRl catalog) 
is used (see §2.3p . 

33. The GALEX NUV effective exposure time in seconds. The non-detections have a 
'-999' entry. For the 991 DR3 quasars covered by a GALEX tile for which no NUV source 
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was found within 7" (a '-999' entry for columns 22-24), the NUV effective exposure time for 
the next nearest near-UV source to the DR3 quasar is given in seconds. 

34. The GALEX FUV effective exposure time in seconds. The non-detections have a 
'-999' entry. Of the 3072 FUV detections, 38 are only detected in the FUV (a '-999' entry 
for columns 19-21), and are not analyzed as true matches in this study (see the discussion 
in 0231) . 

35. The GALEX survey flag. The possible values are integers from 1-4 that specify the 
survey from which the UV identification of the SDSS quasar was drawn (4 = AIS, 3 = NGS, 
2 = MIS, and 1 = DIS; also see §2.11 for the data/survey descriptions). 



3. ANALYSIS 

3.1. Matching Statistics 

Both intrinsic differences among quasars (e.g., due to reddening by dust within the host 
galaxy) and observational effects (e.g., the Ly a forest entering the FUV band at z ~ 0.1 
and the NUV band at z ~ 0.5) will affect the detection probability in UV observations. The 
top panel of Figure [6] illustrates the redshift distribution of the DR3 sample (gray line), the 
NUV detections (blue dashed-line), the NUV-I-FUV (black solid-line) detections, and the 
non-detections (red line), using redshift bins of width Az = 0.1. In the DR3 sample, the 
local minima seen at redshifts z > 2.5 correspond to where the SDSS star-quasar separation 
is less efficient {z ~ 2.7 and z ~ 3. 4), also where quasa rs cross into the stellar locus in 



multi-dimensional color-color space (IRichards et al.l |2002| ). The bottom panel shows the 



respective fraction of DR3 quasars detected for the same distributions using the same color 
scheme. The NUV+FUV detection fraction begins to noticeably decrease as soon as Ly a 
leaves the FUV band {z ~ 0.4). Note that in the NUV band, we recover nearly all of the 
DR3 quasars until z ~ 1.4 (also where the Ly a forest is extinguishing flux in the entire 
GALEX bandpass), with the detection fraction dropping to ~ 50% hj z ^ 2. The resulting 
UV/optical colors from these matches are discussed in §3.21 



For our statistical tests of the matching results, we analyzed the UV detections and 
non-detections by dividing our dataset into four subsamples: (1) the parent sample of DR3 
quasars, (2) those with an NUV detection, (3) those with an NUV-I-FUV detection, and 
(4) those with no UV detection. The FUV-only detections were not analyzed as a separate 
subsample (see discussion in §2.3p . Figures 1^- [TUl display the distributions of each as functions 
of redshift. Galactic extinction in the u band (A^), gf-magnitude, relative optical color A{g — 
i), as well as NUV effective exposure time (tefr)- The inset of each histogram plot also 
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displays the cumulative distribution functions for each parameter. Table H] summarizes the 
results of performing the two-sided Kolmogrov-Smirnov (KS) test on each subsample relative 
to the parent DR3 quasar sample. The results of the KS tests indicate that brightness [g- 
magnitude) and redshift are the primary explanations for non-detections in the UV. Gas and 
dust along our line of sight preferentially attenuates the flux most strongly in the UV bands, 
which is consistent with the KS test results that show Galactic extinction (Au) is likely 
contributing to the lower number of detections in the FUV band. Although the KS statistics 
would suggest teflf as another explanation for the non-detections, an inspection of Figure [TUl 
reveals the subsample of non-detections appears to be more heavily weighted toward larger 
(teff); which is in the opposite sense of what one would expect for the quasars with no UV 
detection. 



3.2. Color-Color Relations 

Using photometry from the DR3 quasar catalog and the GALEX object catalogs, we 
constructed the color-color diagrams in Figures [TT] and [T2l The top panel of Figure [TT] shows 
the locations of SDSS stars (red points and contours) and the quasar sample (black points 
and contours) in — r vs. u — g color-color space. The bottom panel of Figure [11] shows only 
the 3004 quasars with detections in both UV bands that comprise the clean sample, along 
with the same set of SDSS stars. Nearly all of the optically red quasars with {u — g) > 
are not detected in both UV bands. Note that the colors of the DR3 quasars with UV 
detections are uniformly distributed within the same region that the quasar locus in the top 
panel occupies, indicating no strong optical bias in our detections. 

Using the same color scheme as Figure [TT| Figure W]\ shows the UV/optical colors of the 
quasars detected in both UV bands, along with GALEX staxs from selected DIS fields, which 
we matched to the SDSS database. Since nearly all of the DIS fields covering SDSS quasars 
are centered on the high galactic latitude sky, the density of stars on the sky is relatively 
low (especially in the FUV) and most of the FUV detections are near the flux limit. Since 
far fewer quasars have FUV detections and the scatter in the UV photometry is large, the 
regions of {n — u) vs. (/ — n) color-color space occupied by stars and quasars are not as 
well defined. For the bottom panel of Figure [121 the NUV quasar detections were used to 
construct the u — g vs. n — u colors. 

The density of stars in the NUV for a given location on the sky is higher than the 
corresponding density in the FUV by about an order of magnitude, and thus we were able 
to match more GALEX stars to the SDSS database. The star-quasar separation mu — g vs. 
n — u color-color space indicates that the inclusion of UV colors can aid in the color selection 
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of quasars using UV/optical photometry. In addition, the loci of points are consistent with 
reddening along the line of sight through the intergalactic medium — that is to say the bulk 
of the matches do not have "blue" colors in the UV while having "red" colors for the optical 
bands. 



3.3. Color-Redshift Relations 

The observed structure of the quasar colors with redshift is quite complex — in addition 
to intergalactic reddening of the continuum, there are strong emission lines (e.g., Balmer 
recombination lines, Mg II, CIV, and Ly a), and there is the SBB (~ 2000 - 4000 A) due 
mostly to Balmer continuum and Fe II complex emission. Figure [13] shows the color- redshift 
relations for the UV/optical colors f — n, n — u, and g — i colors (all magnitudes corrected 
for Galactic extinction). After binning the quasars by redshift, we have used a cubic spline 
fit to interpolate and subtract the median color at each redshift to obtain the relative colors. 
Broad emission and / or absorption lines will affect the quasar colors as they pass through the 
photometric bands, and are primarily responsible for the observed color-redshift relations. 
For example, in the top-left panel of Figure [13] the Ly a emission line has entered the NUV 
band by z ~ 0.4 and exits at z ~ 1.3, after which the Ly a forest continues to redden 
the UV/optical n — u color with increasing redshift. Quasars that display broad absorption 
lines in their spectra (BALQSOs) can have colors that are either redder or bluer, depending 
on the location of ab sorption troughs relative to the photometric filter bandpasses, and 



Reichard et al.l (12003! ) showed that the overall flux deficit due to reddening dominates over 



the extinction due to the absorption troughs themselves. 

We select the distributions of the relative f—n and relative g—i colors (hereafter A(/— n) 
and A(5f — ^), respectively) shown in Figure [H] for further analysis. T he solid curves a re the 



leas t- Gaussian fits to the blue side and peak of each histogram. [Richards et al.l (120031 ) 



and [Hopkins et al.l (|2004[ ) have shown that the red tail on the A{g — i) distribution is not 
produced by errors in the SDSS photometry, but is rather a true indicator of reddening by 
dust along our fine of sight at the redshift of the quasar. The reddest quasars may not be 
detected in the NUV or FUV, and therefore would not contribute to the observed relative 
color distribution in Figure UM However, the A{g — i) KS-test (^ seems to show that there 
is no significant bias with respect to A{g — i). 

The dashed curve on the top panel of Figure [Tl[ is the result from Monte Carlo simula- 
tions (10^ runs) of A (/ — n) using a Gaussian distribution of photometric errors that use 
the empirical UV magnitude-rms relations as determined in §2.4[ Each best-fit Gaussian has 
a ~ 0.4, which is more than a factor of 2 larger than if the errors quoted in the GRl tables 
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are used. The blue and red curves in the top panel of Figure [T3] indicate the dispersion in 
A {f — n) is dominated by photometric uncertainties, and that it is not purely the result of 
intrinsic differences in the UV quasar SEDs. However, the red tail of the A(/ — n) distri- 
bution (a ~ 0.13) suggests a marginally detected population of quasars heavily reddened in 
the UV, which is also supported by a red tail observed in the A{n — u) distribution. 

In terms of (^-magnitude, Au, and redshift, the quasars with both red A(/— n) and A{g — 
i) colors are statistically indistinguishable from those with bluer relative colors (see Figure 
[TBil . From computing the Spearman Rank-Order correlation coefficients, the significance level 
(~ 0.88) indicates that it is unlikely the relative colors are correlated (or anticorrelated) in 
any way. The quasars with red UV colors in our detections could be the result of intrinsically 
redder SEDs, the Ly a forest, or due to the presence of emission or absorption lines in one 
of the bands that can produce quasar colors that are redder than the bulk of the quasars at 
that redshift. Therefore, while the use of UV colors may aid in the selection o f dust-reddened 



quasa rs, it will also select normal (but perhaps optically steep) quasars (see [Richards et al. 



20031 . and references therein). 



3.4. SEDs 

After applying Galactic extinction corrections to the UV/optical photometry and de- 
redshifting the effective frequency of each filter, we constructed the rest-frame SEDs shown in 
Figure [161 The bluest points are the FUV photometry, with the reddest points representing 
the reddest optical band (the SDSS-z filter). The black line is the median smoothed SED 
for all of the detections, and the vertical dashed lines show the effective frequencies of each 
photometric band, after de-redshifting using the median redshift of all of the detections 
{z = 1.01). Figure [T71 shows the result after each individual SED was normalized to its 
own flux at 2200 A by linear interpolation in log (z^/i/) vs. log (u) space. A cubic spline 
interpolated median SED is the curve shown below the colored points, shifted downward 
1 dex in logarithmic relative flux for clarity. The signatures o f various canonical f eature s 



known to be present in quasar composite spectra (such as that of IVanden Berk et al.l (120011 )) 
are visible, denoted by reference lines given for the locations of Ly a, C IV A1550, the 2800 
A bump (Fe II complex, Mg II, and Balmer continuum emission), H f3, and O [III] A5007. 

In addition to observational effects (e.g., line-of-sight absorption), measurement uncer- 
tainties, and any intrinsic differences in the quasar emission, variability also plays a role in 
the variety of SEDs present in Figure [T71 The ant icorrelations between the variability am- 



plitude and wavelength, as well as luminosity (see IVanden Berk et al.l l2004j . and references 



therein) render UV observations more sensitive to changes in emergent flux levels. Vari- 
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ability can therefore not only displace quasars from their true locations in multidimensional 
color-color space, but also increases the dispersion of the SEDs in our sample. In principle, 
foreground absorption systems can exclude a significant fraction of the quasar population 
from flux-limited sur v eys su ch as the SDSS due to attenuation in the res t-frame UV. How - 



ever, 



Akerman et al.l (120051 ) found no direct evidence for this effect (c.f. lYork et al.ll2006l ). 
The scatter in Figure [17] extending toward larger relative flux at longer wavelengths is pri- 
marily due to the increasing contribution of the host-galaxy to the optical flux for the lower 
redshift quasars, which dominate the SED near the long wavelength limits of our composite 
SED. The peak flux of the SED occurs near ~ 1000 A or so, before the Ly a forest begins 
to significantly extinguish the detected UV flux (Ly a has been redshifted beyond the NUV 
band by z ~ 1.33). The rest-frame wavelength coverage of the UV/optical SEDs extends 
shortward to ~ 350 A for the highest redshift detection at 2; = 3.41 and longward to ~ 9000 
A for the lowest redshift detection at z = 0.12. Table O illustrates the format and content 
of our overall SED that accompanies the electronic edition of the journal. 

In order to constrain our understanding of the UV spectral slope, as well as its red- 
shift and luminosity dependence, one must take into account various observational effects. 
A significant fraction of the quasar luminosity comes from broad emission lines, and Ly- 
man absorption along our line of sight can significantly alter the quasar colors (especially 
beyond z ~ 1). In addition, there is an artificial correlation between quasar redshift and 
luminosity (e.g.. Figure [T8|) . and Lyman absorption strongly attenuates the SEDs. The anal- 
ysis of any luminosity-dependence of the SEDs requires the sample to have a large range 
in luminosity (and therefore redshift), and small redshift windows allow for study of the 
luminosity-dependence as a function of wavelength. In order to remove the effects of the 
emission lines, Lyman absorption and to account for the L — z correlation, we divide our 
UV detections with redshifts 0.11 < z < 3.41 into three redshift bins: (1) 0.15 < log (1 -|- z) 
< 0.20, (2) 0.30 < log {l + z) < 0.35, and (3) 0.45 < log (1 + ^) < 0.50 (see the highlighted 
regions of the L — z distribution in Figure [T5I) . We then computed median SEDs by binning 
the quasars in each redshift bin by their 2200 A luminosity. The resulting median UV/optical 
SEDs are shown in Figure [191 The three panels each correspond to the redshift range in- 
dicated; the blue, red, and green points and lines represent different luminosity bins — the 
median monochromatic luminosity at 2200 A for the quasars in each bin (denoted by Lmed) 
is given in units of log ('^-^2200A'' error bars on each median SED are the la rms 

(the standard deviation normalized by a/ZV, where N is the number of objects) calculated 
for all objects in each bin (gray points). For several points, the errors in the median are 
smaller than the plotting symbols. 

The weak anticorrelation between the UV power-law spectral slope [ay) and the quasar 
luminosity for the lowest redshift bin (bottom panel of Figure [T9l) is primarily due to con- 
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tamination by the host-galaxy, as the lower luminosity quasars are preferentially at lower 
redshifts where the host-galaxy contribution can be comparable to the nuclear emission. A 
discussion of the relations hip between the AGN and host-galaxy fluxes of SDSS broad-line 
active galaxies is given by IVanden Berk et al.l (120061 ) . Our results also reveal a significant 
anticorrelation between a,^ and luminosity for our highest redshift bin (top panel of Figure 
[T9l) covering the shortest rest-frame wavelengths, in the sense that the higher luminosity 
quasars have a steeper UV SED. In terms of the relative FUV flux for each median SED 
in the highest redshift bin, the lowest luminosity {Lmed = 46.5) UV continuum level is 37% 
higher {5.2a) than for the most luminous {Lmed = 47.0) median SED. 

The observed luminosity-dependence of the UV SED sha pe at the short est wavelengths 



is consistent with idea that the Baldwin Effect (hereafter BE; iBa 



dwin IQTtI) mav be du e to 



intrinsic differences in the strength of the ionizing continuum (see lOsmer fc Shields! 1 19991 . for 
a review). Other studies have also drawn similar conclusions, while also taking into account 
the importance of metallicity in dete rmining the emission line strengths as functions of the 
SED shape (e.g.. iKorista et al.lll998l ). For the higher luminosity quasars (having apparently 
redder SEDs), fewer incident UV photons are available to photoionize the broad-line gas 
and will result in lower emission line equivalent widths. Except for the strong luminosity- 
dependence of our SEDs at the shortest wavelengths, and the host-galaxy contributions seen 
at the longest wavelengths, the median SEDs within each redshift bin are largely indistin- 
guishable. 



4. DISCUSSION AND CONCLUSIONS 

We have presented the results of GALEX GRl UV matches to optically-selected quasars 
in the SDSS DR3 quasar catalog. In the final clean sample, 81.1% of the DR3 quasars have 
an NUV detection, while 47.2% have a detection in both the NUV and FUV bands. An 
analysis of the subsamples of the detections and non-detections indicates that the non- 
detections tend to be optically faint objects or at high redshift. The accuracy of the SDSS 
astrometry has allowed us to obtain a rigorous estimate of the accuracy of GALEX positions 
— showing them to be accurate to 0.6-0.7". We have also carried out an empirical analysis 
of the GALEX UV photometry using multi-epoch tiles from the Deep Imaging and All-Sky 
Imaging Surveys that shows the true photometric errors are a factor of ~ 2 larger than those 
given in the GRl object tables. The distributions of UV/optical colors of the quasar matches 
indicate that the inclusion of UV/optical colors can aid in the color selection of quasars, but 
also that large uncertainties in the GALEX UV photometry prevents the use of UV colors 
alone to separate stars from quasars. In addition to other observational effects (e.g., Lyman 
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absorption), extinction within our own Galaxy is likely causing some of the optically-selected 
quasars to be missed by GALEX, especially in the FUV band. Although the observed width 
of the relative UV f — n color distribution is dominated by photometric error, the red tail 
suggests a marginally detected population of reddened quasars — either due to intervening 
absorption systems or due to dust at the quasar redshift. 

In order to investigate a possible luminosity-dependence of the shape of the UV quasar 
continuum, we have constructed rest-frame SEDs for 3 004 quasars — for which we can also 



make comparisons to theoretical models. For example, iHubeny et al.l (120001 ) have used non- 
LTE accretion disk models to predict the emergent SEDs, and their results predict both a 
change in the slope of the power-law continuum in the vicinity of Ly a and a steepening of 
the continuum with increasing quasar luminosity (for fixed values of L/L-Edd)- To properly 
account for the correlation between redshift and luminosity, we have analyzed our SEDs in 
independent redshift and luminosity bins. Our results show an anticorrelation between a^, 
and the quasar luminosity at the shortest wavelengths at the 5.2a level, which appears to 
be consistent with the idea that the BE may be d ue to intrinsic differences in the strength 
of the ionizing continuum. The IScott et al.l (l2004l ) resu lts show a harder U V c ontinuum in 



their F L/^E" composite relative to the if5T composites of IZheng et al. 



Jl997h and 



Scott et all fl2004f ) note 



elfer et al. 



(I2OO2I ). which display a break in the spectral index near ~ 1000 A. 
that their sample of A GN has a median l uminosity that is nearly an order of magnitude 
lower than the sample of iTelfer et al.l (120021 ). The anticorrelation we find in the UV between 
and luminosity is qualitatively consistent with the idea that the Ff/iS'i? composite having 
a significantly harder spectral index (a = 0.56) than the i/ST composite (a = 1.76), is due 
to the lower median luminosity of the FUSE sample. 

Despite the range of L/L^dd spanned by our s ample, our fin d ings a re also in good 
quahtative agreement with the model predictions of iHubeny et al.l (120001). Our observed 
luminosity-dependent UV SED shape is also consistent with the results of lShang et al.l (120051 ) 
making use of both FUSE and HST spectra, who re port a softer FUV s pectral index for the 
higher luminosity AGNs in their sample. However, IShang et al.l (120051 ) also note that most 
of their objects display a spectral break near ~ 11 00 A, which is in rough agreement with 
pre yious work u s ing H ST composite spectra (e.g.. iTelfer et al.ll2002l ). Therefore, although 



the IShang et a 



J2005h result s seem to confirm the change in the continuum slope near Ly 



a predicted by IHubeny et al.l (l2000l ). it remains unclear as to whether or not the location of 
any UV spectral break is a function of luminosity. 

Some work has suggested that the shape of the ionizing continuum may be correlated 
with metallicity, or that it may be the result of luminosity-dependent spectral variations 
(e.g. iDietrich et al.ll2003l ). Differences in the Lyman absorption along different sightlines 
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are unlikely to cause the luminosity-dependence we see, since the composite SEDs are av- 
eraged over a large number of objects at nearly equal redshifts. O ur study has not taken 



into account the presence of intervening absorption systems (e.g., lYork et al.ll2006l ). The 
GALEX photometric bands are wide, covering a significant range of Lyman absorption — 
a rigorous treatment of corrections that take into account absorption by the Ly a forest, 
Lyman limit systems (LLSs), and damped Ly a absorbers (DLAs) is beyond the scope of 
this paper. However, our large sample size has permitted detailed study of the UV/optical 
regime without being highly sensitive to intervening galaxy absorbers or lines of sight that 
intersect LLSs and/or DLAs. Future SDSS and GALEX data releases will permit the con- 
struction of even larger UV/optical samples of quasars. In addition, the larger data volume 
and the anticipated improvements in the photometric calibration of GALEX data released 
in GR2 should increase the number of UV detections of SDSS quasars, as well as aid the 
future studies that make use of UV/optical colors in the color selection of quasars. 



The GALEX early release data was made available in 2003 December. The first major 
data release occurred on 2004 October 1, coincident with the start of the GALEX guest 
investigator program. All GALEX data are served by the Multimission Archive at the Space 
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The grism, imaging window, and uncoated aspheric corrector were supplied by France. We 
acknowledge the dedicated team of engineers, technicians, and administrative staff from the 
Jet Propulsion Laboratory/Caltech, the Orbital Sciences Corporation, the University of Cal- 
ifornia at Berkeley, the Laboratory Astrophysique Marseille, and the other institutions who 
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Fig. 1. — Positional offsets for GALEX GRl matches in both UV bands (3034 quasars) to 
SDSS DR3 quasars. The circles correspond to search radii 9 of CZ/Q, l'/5, 1'.'7, and 2'.'6 that 
recover 68%, 90%, 95%, and 99% of the matches, respectively. The x's identify the 30 objects 
found for 2'/6 < < 7'.'0 that represent the outer 1% of the detections. The small points 
(3004 objects) represent the clean sample and are the quasars referred to and analyzed as 
'detections' throughout the remainder of this paper. The median coordinate offsets {GALEX 
- SDSS) for all of the matches in both bands for right ascension and declination are -0'.'03 
and -0'.'12, respectively. 
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Fig. 2. — A comparison of multi-epoch photometric measurements for the GALEX Deep 
Imaging Survey (DIS). A total of 6965 objects remain after outliers > Sa from the median 
value for the n-magnitude differences have been removed — likely to be spurious matches 
between both epochs of the same GALEX tile. The red lines are the photometric errors 
taken from the GRl merged object catalogs for each object after adding the errors from each 
epoch in quadrature, and are reflected about An = 0; the blue curve represents the intrinsic 
rms dispersion calculated using appropriately sized n-magnitude bins and is by definition 
symmetric about the n-magnitude axis; the green curves are the calculated 68.3% confidence 
limits for the objects. The intrinsic rms very nearly overlaps our 68.3% confidence limits — 
indicating that the photometric errors are approximately Gaussian. The figure suggests that 
the true photometric errors are roughly a factor of ~ 2 larger than the published ones. 



- 25 - 




Fig. 3. — Using the same notation as Figure[2]for the GALEX All-Sky Imaging Survey (AIS) 
observed at multiple epochs. 3966 objects remain after outliers > ?>a from the median value 
for the n-magnitude differences were removed. 
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Fig. 4. — The binned rms magnitude errors as a function of magnitude from matches on 
multi-epoch GALEX tiles for the NUV (top) and the FUV {bottom). The sohd and dashed 
lines represent the least-x^ fits to the binned distribution of empirically determined photo- 
metric errors, for AIS and DIS fields respectively. 
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Fig. 5. — Histogram showing the overall distribution of NUV effective exposure times for 
the DR3 quasars detected by GALEX. For the objects to the left of the empirical division at 
log (tcff) =2.7 {dashed line), we use the AIS empirical curves of Figure HI and for objects to 
the right of this line, we use the DIS empirical curves to compute our estimates of the true 
NUV and FUV photometric errors for all of the UV detections in our sample. 
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Fig. 6. — Quasar redshift distributions, (top) — The sohd black hue shows the binned 
distribution of quasar redshifts (0.12 < z < 3.41; median 1.02) for the 3004 G ALEX /SDSS 
detections in both UV bands. The same distribution for the DR3 quasars covered by all 
GALEX GRl fields (0.12 < z < 5.20; median 1.45) is the solid gray line. The blue dashed 
line is for all of the DR3 quasars with an NUV detection, and the red line represents quasars 
without a UV detection. The local minima in the DR3 distribution (at z ~ 2.7 and at 
z ~ 3.4) correspon d to redshifts wh e re the DR3 target selection algorit hm has lower efficiency, 
as is discussed by ISchneider et al.l ( 120051 ) and [Richards et al.l (120061 ) . The inset shows the 
cumulative distribution functions for each subsample, using the same color scheme, (bottom) 
— The detection fraction for the quasars found in both GALEX bands (black solid line) as 
well as for the NUV (blue dashed line). In the NUV, nearly all of the DR3 sample quasars 
out to 2; ~ 1.7 are detected by GALEX. 
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Fig. 7. — Using the same notation as Figure El the distributions of Galactic extinction in 
the SDSS M-band using bins of width AAu = 0.025. 
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Fig. 8. — Using the same notation as Figure [6l here for the optical (^-magnitude distributions 
(no Galactic extinction corrections have been applied here — see using magnitude bins 
of width Ag = 0.1. 
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A(g-i) 

Fig. 9. — Using the same notation as FigureEj the relative color A{g — i), which is defined as 
the difference between the g — i color of the quasar and the median g — i color at its redshift. 
The bins have a width of 0.1. 
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Fig. 10. — Using the same notation as Figure El here for the NUV effective exposure time 
(tcfr) as defined in §2.31 using bins of width dlog tefr = 0.2. 
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Fig. 11. — {top) The black x's and contours show the distribution of colors for 6371 SDSS 
DR3 quasars covered by GALEX GRl fields. The red points and contours are stars from 
the DR4 database, (bottom) Identical to the top panel, but only the quasars detected in 
both GALEX UV bands are shown along with the same set of SDSS stars. All objects have 
been corrected for Galactic extinction. 
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Fig. 12. — The black points and contours show the distribution of colors for 3004 quasars 
detected in both GALEX bands. The red points are stars from GRl DIS fields matched to 
DR4 stars {g < 20.5). All objects have been corrected for Galactic extinction. 
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Fig. 13. — UV/optical colors as a function of redshift for all of the UV detections. The 
median color in each redshift bin was interpolated for each quasar at each redshift (black 
curves). Galactic extinction corrections have been applied to all objects. 
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Fig. 14. — The UV/optical relative colors for all of the UV detections, computed by sub- 
tracting the median color at each redshift. The solid curves are the least-squares best-fit 
Gaussians to the peak and blue tail of each histogram. Note the red tails of both distribu- 
tions that rise above the fits (see §3.30 . The red tail in l\{g — i) is not produced by the SDSS 
photometric uncertainties, but is rather due to reddening by dust alon g our line of sight or at 
the redshift of the quasar (IRichards et al.l l2003l : iHopkins et al.l |200J) . The dashed curve in 
the top panel is the result of Monte Carlo simulations of the A(/ — n) distribution using our 
empirical UV photometric errors ( §2.4p . The comparable widths of both Gaussians indicates 
that the broader distribution of / — n colors is not due to intrinsic differences in the UV 
SEDs of the detections, but is dominated by the photometic uncertainties. However, the 
'excess' population of objects with large A(/ — n), may suggest the presence of a subset of 
heavily reddened objects. 
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Fig. 15. — The points and contours show the UV/optical relative colors for the quasar 
detections on both GALEX bands. The Spearman Rank-Order Correlation results do not 
indicate any significant correlation (or anticorrelation) of the UV/optical relative colors (see 
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Fig. 16. — Rest-frame spectral energy distributions (SEDs) for all of the detections (see 
§2.3p . There is a one-to-one correspondence between the color of the plotting symbols and 
the effective wavelength of each band (e.g., the FUV is violet and the z-band is red). The 
solid curve is the smoothed median SED using dlog(z/) = 0.01. The dashed gray vertical 
lines correspond to the indicated photometric bands, shown with the effective frequency of 
each de-redshifted using the median redshift {z ~ 1.02) for all of the objects. 
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Fig. 17. — Our rest-frame SEDs from Figure [16] where each individual SED has been nor- 
mahzed to its own flux at 2200 A. The peak is near the location of Ly a (1216 A), where 
the overall rise in the relative flux is largely due to the presence of the Ly a emission line 
within the filter bandpasses. The locations of other prominent features are also visible, such 
as C IV 1550, the 2800 A bump, as well as H/5 and O [III]. 
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Fig. 18. — The specific luminosity (at 2200 A; rest-frame) vs. redshift distribution of the 
3004 NUV/FUV GALEX detections of DR3 quasars. The three highlighted redshift ranges 
correspond to 0.15 < log [l + z) < 0.20, 0.30 < log (l + z) < 0.35, and 0.45 < log (1 + ^) < 
0.50, respectively. We further divide the quasars within each redshift bin into three bins of 
luminosity, which are centered on the median value and with each luminosity bin containing 
a third of the quasars. We have used the quasars in each bin to construct the UV/optical 
SEDs of Figure M 
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Fig. 19. — SEDs in bins of redshift and luminosity. The color points and lines are the median 
normalized SEDs (see Figure [T7j) constructed from independent redshift and luminosity bins, 
drawn from the luminosity-redshift distribution of Figure [18] over the indicated redshift 
range. The median log (z/Lj^) at 2200 A of the quasars used to construct each individual 
SED are given by Lmed- In the top and middle panels, Lyman absorption and extragalactic 
extinction are the primary reasons for the SED turnover near log u ~ 15.4 (Ly a). The error 
bars (comparable in size to the plotting symbols) are the la rms (the standard deviation 
normalized by \/iV, where N is the number of objects in each bin) computed for all of the 
quasars (gray points) contributing to each respective SED. The host-galaxy contamination 
becomes increasingly important at longer wavelengths, especially at lower redshifts where 
the correlation is easily visible in the bottom panel. In the top panel, our SEDs cover the 
shortest rest-frame wavelengths, where we see a strong anti-correlation of the overall strength 
of the UV continuum with luminosity. 
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Table 1. Best-Fit Parameters for Empirically Determined Photometric Errors 



Survey Band a b 

AIS NUV 0.20525 -4.7302 

AIS FUV 0.13360 -3.0578 

DIS NUV 0.23466 -5.7196 

DIS FUV 0.24413 -5.8804 



Note. — Parameters for the least- 
fits to the binned NUV and FUV 
rms as a function of NUV and FUV 
magnitude using Equation [1] in Figure 
il 
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Table 2. Table for GALEX GRl Matches to SDSS DR3 Quasars 



Object (SDSS J) 


GALEX Object I.D. 


R.A. (SDSS) 


Decl. (SDSS) 


Coord. Offset (") 


Redshift 


000009.42-102751.9 


2435007175665396402 


0.039269 


-10.464428 


1.62 


1.8442 


000011.96+000225.3 


2667127184236220309 


0.049842 


0.040372 


0.44 


0.4790 


000017.38-085123.7 


2434901622551226166 


0.072421 


-8.856607 


1.17 


1.2499 


000026.29+134604.6 


2665121679322122922 


0.109578 


13.767970 


0.65 


0.7678 


000028.82-102755.7 


2435007175666450605 


0.120086 


-10.465496 


0.32 


1.1377 



Note. — A full version of the above table accompanies the electronic edition of the journal. Table [2] serves 
as a reference for its hnal form and content. All object coordinates are given in decimal degrees (J2000.0). 
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Table 3. Table Format for GALEX GRl Matches to SDSS DR3 Quasars 



Column Data Format Descrii:)liou 



1 


A18 


bubb Diio Quasar Catalog object designation m nnmmss.ss+ddmmss.s (JzOOO.Oj 


2 


A19 


GALEX GRl merged object catalog designation; -999 indicates no detection 


3 


b 11.6 


bDSS right ascension m degress (J2000.0) 


4 


F11.6 


SDSS dechnation in degrees (J2000.0) 


5 


F11.6 


GALFa right ascension (degrees) 


6 


F11.6 


GALEX dechnation (degrees) 


7 


F7.2 


Coordinate onset between GALEX/SDSS positions given m arcseconds 


8 


F7.4 


redshift (SDSS) 


9 


F7.3 


SDSS u magnitude (not corrected lor Galactic extinction) 


10 


F6.3 


error m hDSb Fbl< u magnitude 


11 


F7.3 


SDSS g magnitude (not corrected tor Galactic extinction) 


12 


F6.3 


error m SDSS FSl" g magnitude 


13 


F7.3 


SDSS r magnitude (not corrected lor Galactic extinction) 


14 


F6.3 


error m SDSS PSF r magnitude 


15 


F7.3 


SDSS I magnitude (not corrected tor Galactic extinction) 


16 


F6.3 


error in SDSS PSF i magnitude 


17 


F7.3 


SDSS z magnitude (not corrected lor Galactic extinction) 


18 


F6.3 


error m SDSS rbb z magnitude 


19 


F9.3 


GALEX / magnitude (not corrected lor Galactic extinction); -999 ii no detection 


20 


F8.3 


error in PSF / magnitude (from GRl object catalog); -999 if no detection 


21 


F8.3 


error m FSl^ / magnitude (empirical calibration) 


22 


F9.3 


GALEX n magnitude (not corrected for Galactic extinction); -999 if no detection 


23 


F8.3 


error in PSF n magnitude (from GRl object catalog); -999 if no detection 


24 


F8.3 


error in PSF n magnitude (empirical calibration) 


25 


F7.3 


Galactic extinction in the u band {Ay) 


26 


F7.3 


Logarithmic Galactic H I column density 


27 


F8.3 


Mi {Ho = 70 km s"^ Mpc'^f^M = 0.3, = 0.7,0;^ = -0.5)* 


28 


F8.3 


Ljy at 2200A (rcst-framc; -999 if no coverage) 


29 


F8.3 


Ljy at 5100A (rcst-framc; -999 if no coverage) 


30 


F6.3 


A{g — i) relative color (corrected for Galactic extinction) 


31 


F8.3 


A(/ — n) relative color (corrected for Galactic extinction) 


32 


F7.5 


object distance from GALEX held center (degrees) 


33 


F8.3 


GALEX near-UV effective exposure time (seconds) 
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Table 3 — Continued 



Column 


Data Format 


Description 




34 


F8.3 


GALEX far-UV effective exposure time (seconds) 




35 


11 


GALEX survey identifier (4 = AIS, 3 = NGS, 2 = 


: MIS, 1 = DIS) 


Note. 


— Data format 


and content for the full version of Table [2] that 


accompanies the 



electronic edition of the journal. 



Table 4. Kolmogorov-Smirnov Test Results Summary 



Sample 1 


Sample 2 




(^-magnitude 


A(<7-i) 


Redshift 


NUV teff 


NUV+FUV detections 


DR3 Sample 


0.18 


< 10-12 


0.69 


< 10-12 


< 10-12 


NUV detections 


DR3 Sample 


0.80 


4.3 X 10-12 


0.67 


< 10-12 


< 10-12 


no UV detection 


DR3 Sample 


3.8 X lO"-'^ 


< 10-12 


0.63 


< 10-12 


< 10-12 


NUV+FUV detections 


no UV detection 


2.5 X IQ-^ 


< 10-12 


0.27 


< 10-12 


< 10-12 


NUV detections 


no UV detection 


3.4 X 10"^ 


< 10-12 


0.62 


< 10-12 


< 10-12 



Note. — Our Kolmogorov-Smirnov significance parameters computed for the indicated subsample dis- 
tributions of selected parameters, after matches with 9 > 2'.'6 have been removed. 
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Table 5. Sample Table for Median SED 



Log A (A) 


Log V (Hz) 




rms 


3.831 


14.646 


0.810 


0.041 


3.788 


14.689 


0.699 


0.030 


3.752 


14.725 


0.611 


0.018 


3.731 


14.747 


0.610 


0.026 



Note. — A full version of the above 
table accompanies the electronic edition of 
the journal. Table [5] serves as a reference 
for its final form and content. 



